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We present numerical results for an S = 1/2 Heisenberg antiferromagnet on a inhomogeneous 
square lattice with tunable interaction between spins belonging to different plaquettes. Employing 
Quantum Monte Carlo, we significantly improve on previous results for the the critical point sepa- 
rating singlet-disordered and Neel-ordered phases, and obtain an estimate for the critical exponent u 
consistent with the three-dimensional classical Heisenberg universality class. Additionally, we show 
that a fairly accurate result for the critical point can be obtained from a Contractor Renormalization 
(CORE) expansion by applying a surprisingly simple analysis to the effective Hamiltonian. 
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I. INTRODUCTION 

The square lattice quantum Heisenberg antiferromag- 
net (SLQHA) with spin S" = 1/2 is one of the paradig- 
matic models in condensed matter physics and has been 
extensively investigated in the last two decades, mainly 
in connection with the parent compounds of cuprate 
supcrconductorsji Good agreement with experimental 
data is obtained from the analysis of an effective con- 
tinuous field theory, rigorously justified in the limit of 
large spin S, given by the (2 + l)-dimcnsional nonlin- 
ear a model (NLcrM).— The coupling g in the NLcrM 
controls the transition between Neel-ordered and quan- 
tum disordered phases and it can be shown that the 
SLQHA maps to the renormalized classical regime of the 
NLcrM {g < gc), which therefore has a long-range ordered 
ground-statcjii^ 

The proximity to a quantum critical point may, how- 
ever, lead to a crossover to a finite temperature quantum 
critical regime where thermodynamic quantities are af- 
fected by strong quantum fluctuations and display uni- 
versal behaviori^ Nevertheless, numerical simulations^ 
indicate that this crossover may be too narrow to be de- 
tectable, possibly due to the fact that the SLQHA is deep 
inside the renormalized classical regime. This difficulty 
has motivated the investigation of antifcrromagnets de- 
fined on decorated square lattices where the proximity 
to a quantum critical point separating singlet-disordered 
and Neel-ordered phases is controlled by adjusting cou- 
plings in the Hamiltoniani^i^iii^i^ii^ 

We consider an S = 1/2 Heisenberg Hamiltonian de- 
fined on the "plaquettized" square lattice depicted in 

Fig.m 
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J {{i,j), bold hues in Fig. [T|) and J' {{i,j) , thin lines 
in Fig. [T]) are, respectively, intra- and inter-plaquctte 
nearest- neighbor antiferromagnetic interactions. Since 
the model is self-dual under the transformation J ^^ J' , 
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FIG. 1: (Color online) The "plaquettized" square lattice con- 
sidered in this paper: nearest-neighbor spins lying on the 
same (neighboring) plaquette(s) interact via superexchange J 
(J'), represented by thick (thin) continuous lines [see Eq. fl])]. 
Dashed lines highlight the clusters employed in obtaining the 
range-1 (a), -2^''^ (b) and -2 (c) CORE results. 



we apply the restriction J' < J without loosing generality 
and set J = 1. Similarly to what happens with the afore- 
mentioned models with tunable interactions^^iiiSiiSiiiS 
the ratio J' /J (equivalent to g~^ in the NLcrM) con- 
trols the magnitude of quantum fluctuations and a quan- 
tum phase transition at Jq/J, belonging to the uni- 
versality class of the 3D classical Heisenberg model)^ 
separates a disordered singlet phase at low J'/ J from 
the renormalized classical state at J' /J = 1, where 
the original SLQHA is recovered. Previous results for 
this "plaquettized" Heisenberg antiferromagnet were ob- 
tained analytically^ from series expansions^^^ exact 
diagonalization of small clustersHiH and by diagonalizing 
the effective model obtained from a Contractor Renor- 
malization (CORE) expansion F^^ii^ The currently best 
estimate for the critical point [Jq/J ~ 0.555(10)] was 
obtained from an Ising series expansionJ^ 

We investigate the model defined by Eq. ^ by means 
of Quantum Monte Carlo (QMC) simulations and CORE 
with a twofold purpose. Firstly, we would like to im- 
prove on previous estimate o-'^-'^'-'^^i-'^^'-'^^i"'^^'"'^^ for the quan- 
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FIG. 2: (Color online) (a) SSE-QMC data for spin stiffness 
ps multiplied by the system's size L as a function of the inter- 
plaquette coupling J' / J, for temperatures T < 1/21/. Error 
bars are much smaller than the symbols' size, (b) Conver- 
gence of the intersection points for curves L and 2L for psL 
[data shown in (a)] and second-order Binder cumulant Q2 (not 
shown). Extrapolation (see main text) gives us the estimate 
J'c/ J ~ 0.5485(12) for the critical point, (c) Data collapse for 
psL is attained for Jq/J = 0.5472(8) and u = 0.71(1), using 
the scaling Ansatz of Eq. (Q. 



turn critical point. This might pave the way to future 
investigations of, for instance, the effects of impurities 
in a quantum critical antiferromagnet.-^^^'^ Secondly, 
we are interested in testing the quality of the results ob- 
tained from a CORE expansion including longer-ranged 
effective interactions than in Refs. llSlfla and by applying 
a simpler analysis to the effective Hamiltonian. 



II. NUMERICAL RESULTS 
A. Quantum Monte Carlo 

We have performed QMC simulations for the model 
described by Eq. ([1]) by employing the ALPS libraries^S 
implementation of the directed loops algorithm^ii^ for 
the Stochastic Series Expansion (SSE) representation.— 
Lattices with L x L sites {L/2 x L/2 plaquettes) with L 
up to 36 have been considered, with periodic boundary 
conditions along both directions. Temperatures are set 
so to ensure that ground-state properties are accessed: 
for each system's size, simulations were performed for 
inverse temperatures /3„ — 2" so that the obtained av- 



erages agreed (within error bars) for /?„ and /3„ 
have calculated the spin stiffness ps and the second-order 
Binder cumulant for the staggered magnetization, Q2- 
The spin stiffness is obtained in terms of the winding 
numbers w^ and Wy, 
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and is expected to scale close to the critical point like 
Ps ^ L'^~'^~^, where d = 2 is the dimensionality and 
the dynamic critical exponent is expected to be z = 
1,^ Therefore the quantity psL should assume a size- 
independent value at the critical point, something con- 
firmed by our results shown in Fig. [H^a). The second- 
order Binder cumulant for the staggered magnetization 
(77ig) is defined as 
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Q2 also displays universal behavior in the critical regime 
and curves for different lattice sizes cross close to the 
critical point (not shown). 

In order to locate the quantum critical point, J'c/ J, 
we first analyze the scaling behavior of the intersection 
points between curves for psL and Q2 obtained for lat- 
tice sizes (L, 2L). Crossing points are determined by 
performing linear and quadratic fits to different data sub- 
sets and deviations between different estimates are used 
in setting (generous) error bars. The so obtained results 
arc plotted in Fig. [H^b) as a function of l/£. We re- 
mark, similarly to what was found in Ref. llO for dimer- 
ized magnets, that fastest convergence to the thermody- 
namic limit is attained for p^L and therefore focus on the 
results for this quantity in what foUowsi^ Since the cur- 
vature for the crossing points decreases with increasing L 
[Fig.ISJb)], an upper bound for J^/J is simply obtained 
by applying a linear extrapolation to the three largest 
(L, 2L) crossing points; a lower-bound is directly given 
by the crossing point for the largest pair (L, 2L). In this 
way, we arrive at the result Jq/J = 0.5485(12).— 

In trying to achieve higher accuracy, and additionally 
estimate the critical exponent associated to the correla- 
tion length, v, we employ the scaling Ansatz 



ps{t,L) = L-'f,,{tL'/n, 
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with reduced coupling t = (J' — Jq)/Jq. By plotting 
PsL versus tL^/^ , and adjusting the values of Jq/ J and 
I/, we achieve data collapse for Jq/ J = 0.5472(8) and 
V = 0.71(1), as shown in Fig.HJc),^ 

Our results for Jq/J are consistent with earlier 
estimates , ^^i^^'^'^i^^'^^i^^ but improve on the previously 
best result [J^/J = 0.555(10), Ref. [l^ by one order of 
magnitude. Furthermore, our estimate for the critical 
exponent v is compatible with the most accurate result 
for the 3D classical Hcisenberg model [v = 0.7112(5), 
Ref. [23], as expected from the mapping onto a NLcrMi^ 
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FIG. 3: (Color online) (a) Spin-gap in the disordered phase of 
the spin model on the modulated square lattice described by 
Hamiltonian Eq. ([1} . Results have been obtained from a pla- 
quette series expansion^^ (circles, Ref. Il2 ) and various range 
CORE analysis (see main text). The vertical dashed lines 
indicate values of J' j J consistent with the critical point ob- 
tained from QMC simulations {J'dJ € [0.5466, 0.5497]). The 
zoom in the inset (b) shows that values for the critical point 
consistent with the ones obtained from QMC are obtained 
from range-2 (squares) CORE results (range-2^" results are 
represented by red circles). 



Further improvements in the values for J'q/ J and v may 
be achieved by simulating larger systems and/or applying 
a more sophisticated data analysis, taking into account 
subleading finite-size corrections i^°i^^ 



B. Contractor Renormalization 

The Contractor Renormahzation (CORE) method^i^ 
is a tool in deriving low-energy effective Hamiltonians 
for lattice models and was previously applied to the spin 
Hamiltonian on the modulated square lattice [Eq. ([T])] by 
Capponi et al.i^'^ We extend their results by deriving a 
longer-ranged CORE expansion. We also notice that our 
CORE expansion is essentially a strong coupling (on-site 
repulsion C/ — > cx)) version of the one derived in Ref. |32| . 
where an effective Hamiltonian for the Hubbard model 
on the square lattice was obtained, but we have the ad- 
vantages of a more natural motivation for choosing the 
plaquettes as elementary blocks and of the absence of 
charge degrees of freedom. 

We start by noticing that the original spin model de- 
scribed by Eq. ((T|) is particularly amenable to a CORE 
analysis, for the strongly coupled plaquettes are a natu- 
ral choice as the elementary blocks (see Fig. [TJ for details 
on the CORE procedure the reader is referred to Refs.ilg 
and [32). The retained low- lying block states are the pla- 
quette's singlet ground-state |s) and the triplet states 
i"), with a = —1,0,-1-1 denoting the total S^ compo- 
nent. This choice for the restricted local basis is justified 



by the fact that \s) and |i") arc the lowest block's eigen- 
states and also, as shown by Capponi et ai^^^ by their 
large weight in the density matrix of a plaquette embed- 
ded in a larger cluster. Effective couplings between these 
retained block states are obtained by subsequently di- 
agonalizing a cluster comprised of connected plaquettes: 
a matching number of cluster's low-lying eigenstates are 
projected upon the basis formed by the tensor products 
of local singlets and triplets, an effective Hamiltonian 
being obtained by imposing the constraint that the clus- 
ter's low-energy spectrum is exactly reproduced and by 
subtracting interactions previously obtained from clus- 
ters involving lesser plaquettes. We employ the clusters 
highhghted in Fig.[T]and label the results according to the 
longest range effective couplings obtained at each step of 
the CORE expansion: range-1 interactions arc obtained 
from the cluster comprised by two connected plaquettes 
depicted in Fig. [TJa), range-2^/^ from a four-plaquette 
cluster [Fig. [T]Jb)] and range-2 from the cluster display- 
ing three ahgned plaquettes shown in Fig.mc). 

The effective model resulting from the above proce- 
dure is expected to evidence dominant microscopic mech- 
anisms at play in the original model and physically sound 
results are ideally obtained by means of a simplified sub- 
sequent analysis. We stress that our approach differs 
from the previous on o^^'^^ in a crucial way: while Cap- 
poni et ali^'^ restricted their CORE expansion to the 
shortest-range and studied the resulting range-1 effective 
Hamiltonian by means of exact diagonalizations, we ex- 
pect that good results arc obtainable in a simpler way 
from an extended CORE expansion including longer- 
ranged interactions. Accordingly, we locate the quantum 
critical point by determining the value of J'/ J where the 
gap for triplet excitations (A) vanishes and estimate A 
from the effective CORE Hamiltonian simply as 



fit - Hti + t2 + ts) 



(5) 



fit is the chemical potential for triplet excitations above 
the singlet ground-state for low J' /J and the second term 
accounts for the triplons' kinetic energy: ii is the nearest- 
neighbor (NN), i2 the next-NN and ta the third-NN hop- 
ping amplitudesi^ In other words, Eq. ([5]) is the energy 
of an isolated triplon in a singlet sea. 

Results for A(J'/J) obtained from range-1, -2^/^ and 
-2 CORE expansions are shown in Fig. [31 compared with 
the ones obtained from the plaquette series expansion 
derived in Ref. [l2|i^ As expected, all results mutually 
agree in the limit of small J' / J, where they are essen- 
tially exact; for larger J'/J, range-2^/^ and -2 CORE 
underestimate A. However, the values of J' /J where 
range-2^/^ and -2 results for A vanish are in surpris- 
ingly good agreement with the QMC results for J^/J 
presented in Sec. IIIAI J'/ J « 0.5513 for range-2i/2 and 
J' /J w 0.5491 for range-2. These values are seemingly 
converging very fast to a result consistent with the QMC 
estimates and this suggests that the procedure employed 
here might lead to more precise estimates for the critical 
point than the one employed in Refs. |l5| and |l6l , where 



J'q/ J ~ 0.55(5) was obtained. We conjecture that this 
unexpected accuracy is related to level-crossings observed 
close to the point where A vanishes. However, obviously 
we cannot discard the possibility that this remarkable 
agreement is coincidental and remark that poorer results 
are obtained for intermediate values of J' / J ^ with the 
consequence that A from CORE does not follow a power 
law (as seen from a logarithmic plot, not shown). How- 
ever, it is desirable to further test the procedure employed 
here by considering similar spin models f 



III. CONCLUSIONS 

Summarizing, we have investigated an 5 = 1/2 Heisen- 
berg antifcrromagnet on a "plaqucttized" square lattice 
(Fig. [H) by means of QMC and CORE. Our results 
for the quantum critical point separating the gapped- 
singlct and Nccl-ordered phases, J'q/ J = 0.5485(12) and 
Jq/J = 0.5472(8), obtained from QMC simulations, sub- 
stantially improve on previous estimateS j-'^-'^i-'^'^'-'^'^i-'^^'-'^^i^^ 
and the obtained critical exponent v = 0.71(1) is consis- 
tent with the tridimensional classical Hcisenbcrg model 
universality clasSf21 as expected from the mapping to a 
NLcrM,^ 

We also highlight the surprisingly good result for the 



critical point extracted from a simple analysis of range- 
2^/^ and -2 effective CORE Hamiltonians. However, it 
is not presently possible to exclude the possibility that 
the good agreement with QMC is coincidental and it 
would be interesting to further test the procedure em- 
ployed here. The fact that CORE is immune to the infa- 
mous sign problem open interesting research possibilities, 
and fcrmionic systems on a geometry similar to the one 
considered here^ may be investigated. 

Note: While preparing this manuscript, and after fin- 
ishing our simulations, we became aware of work by Wcn- 
zel et al^^^ in which dimerized and quadrumerized two- 
dimensional antiferromagnets are investigated by QMC. 
By employing a more sophisticated data analysis and 
simulating larger lattices, they obtain more precise es- 
timates for the critical point for the model considered in 
the present work. 
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